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Abstract 

Hamilton-Jacobi equation for Brans-Dicke theory is solved by using a long- 
wavelength approximation. We examine the non-linear evolution of the inho- 
mogeneities in the dust fluid case and the cosmological constant case. In 
the case of dust fluid, it turns out that the inhomogeneities of space-time 
grow. In the case of cosmological constant, the inhomogeneities decay, which 
is consistent with the cosmic no hair conjecture. The inhomogeneities of the 
density perturbation and the gravitational constant behave similarly with that 
of space-time. 



1 Introduction 



It is generally believed that the consistent theory of quantum gravity is correctly 
described by the superstring theory. The classical theory of gravity is nothing but a 
low energy effective theory of the superstring theory. The classical theory which is 
predicted by the superstring theory has the form of a scalar-tensor theory, the scalar 
field is the so-called dilaton. Therefore, it is necessary to consider the consequences 
of this extra scalar field at least in the phenomena close to the Planck scale. If 
the dilaton field acquires a large mass due to unknown dynamical mechanism, there 
will be no observable macrospic difference between the superstring predicted theory 
and Einstein's general theory of relativity. However, recently, it is pointed out the 
possiblity of the massless dilaton. [||] If so, it is important to study scalar-tensor 
theory more seriously. The simplest scalar-tensor theory is Brans-Dicke theory 
where the dilaton field acts like a dynamical gravitational constant. 

On the other hand, to circumvent the graceful-exit problem of old inflation, 
Brans-Dicke theory is renewed in the inflational universe scenario as extend infla- 
tion. Moreover, Bellido et.al.0 investigate the stochastic inflation formalism in the 
context of Brans-Dicke theory. They described the inhomogeneous universe with 
fluctuations of the gravitational constant. The main idea of stochastic inflation is 
to solve the equations for the inhomogeneous fields in de-Sitter space by separating 
both the gravitational and scalar fields in short wavelength quantum fluctuations, 
which oscillate on scales smaller than the Hubble radius, and long wavelength fluc- 
tuations which are treated as classical fields. Salopek and Bond[|J have developed 
the formalism to treat the long wavelength fluctuations in terms of Hamilton- Jacobi 
equation, which is applicable not only to the inflational theory but also the late 
stage evolution of the density fluctuations as far as the typical scale of the fluctua- 
tion exceeds the Hubble radius. 

The so-called long wavelength approximation has rather long history dating back 
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to Lifshitz and Khalatonikov. 0] Later, Tomita developed the above approxima- 
tion as the Anti-Newtonian scheme. Recently, Salopek and co-orthors elegantly 
formulated the long wavelength approximation in the context of Hamilton- Jacobi 
theory. 0, |8j The direct method of Comer et. al.|| is also useful to calculate the 
higher order correction. Attempts to apply the formalism to the inflationary theory 
and the higher dimensional theory also exist. [|l(J 

In this paper, we shall apply the long wavelength approximation to Brans-Dicke 
theory. We are interested in the non-linear evolution of the long-wavelength inho- 
mogeneities. As the inhomogeneties of the Brans-Dicke field implies the inhomo- 
geneities of the gravitational constant, it is important for the astrophysical phenom- 
ena. In our modest study, of course, we do not intend to make definite statements to 
the astrophysics. However, it is important to investigate how the inhomogeneities 
of the gravitational coupling constant evolve. We start in Sec. 2 by writing the 
Hamilton- Jacobi equation for Brans-Dicke theory. The spatial gradient expansion 
is performed in Sec. 3. As a matter field, the cosmological constant is interesting. 
This case is related to the cosmological no-hair conjecture. Sec. 4 is devoted to these 
subjects. In Sec. 5, the discussions of the various problems are presented. In the 
Appendix, the results of the direct method are explained. 

2 Hamilton-Jacobi Equation for Brans-Dicke The- 
ory 

For simplicity, we will consider the dust fluid matter. Action for the Brans-Dicke 
theory with the dust matter x is given by 

/ = J d'x^l^R - ^cVR0 - ^-(g^d.xduX + m 2 )} , (1) 

where n is the Lagrange multiplier and m is the particle mass which is normal- 
ized to unity below. Here, u is the parameter of the theory. The Brans-Dicke 
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field 4> is considered as the effective gravitational coupling constant. The Hamilton- 
Jacobi equation for Brans-Dicke theory is obtained using the Arnowitt-Deser-Misner 
(ADM) formalism in which the space-time is foliated by a family of space-like hy- 
persurfaces. In the ADM formalism, the metric is parametrized as 

ds 2 = -N 2 dt 2 + ^{dx 1 + N i dt)(dx j + N j dt) , (2) 

where N and N* are the lapse and shift functions respectively, and 7^ is the 3-metric. 
Using the above metric, we obtain the Hamiltonian form of the action as 

I = j d 4 x{7T ij 'j ij + 7^0 + 7T X X ~NH - WHi} , (3) 

where 

1 a ku UJ + 1 i , 1 ^2 11 



+ V 1 + T j X,iX,j * x - y/l<t>R + u^'diWrf + 2^A0 (4) 
Hi = -2{ 1%k n^),+rc lk llk ,TiU, + ^ x X,- (5) 

Here 7r y , 7r^ and n x are conjugate to 7^,0 and Xi respectively, and R denotes the 
3-dimensional scalar curvature. Variation with respect to the momentum yields the 
equations of motion 

= 2^-3^-^ (6) 
^(X-N 1 X>1 ) = y/l+yix,iXj , (7) 

Variation of the action (3) with respect to the field variables yield the evolution 
equations for the momentum. They are automatically satisfied provided that 

n tJ = -=— , tt* = — , tt x = — 9 

hij 09 ox 
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satisfy the constraint equations and provided the evolution equations (6) hold. Here, 
instead of solving the equations of motion for the momentum fields, we will use the 
Hamilton- Jacobi method. Hamilton- Jacobi equation is 
1 5S 5S _ uj + 1 ^ 

'itikljl ~ : ^zlijlkl] 



1 5S 2 _l 1_ 5£^5S_ 

+ 2(2cj + 3) y/T5<l>' 2u; + 3 ^f lij S-fij Scj) ^ ' 

& s / — 

The momentum constraint is rather trivial condition which states that the generating 
functional is invariant under the spatial coordinate transformation. 

Hamilton- Jacobi formalism has a great advantage, that is, its intimate relation 
to the quantum gravity. Wheeler-DeWitt equation for Brans-Dicke theory is given 
by 

= , = , (11) 

where the canonical commutation relations 



ivix) ,n»(v)\ = l -(6*6 l j + 6W)6(x-y), (12) 
'<j>(x) ,*+(yj\ = iS(x-y) , (13) 
[ X (x) ,nx(y)\ = iS(x-y) (14) 

are assumed. If we consider the WKB approximation, we get eq.(10) as the lowest 
equation. The research in this direction from the point of view of the long- wavelength 
approximation is the issue to which we shall attack in the future. 

3 Long Wavelength Solution 

The heart of the long wavelength approximation is the following: For illustration, 
we take metric in the synchronous form 

ds 2 = -dt 2 + -f ij (x k ,t)dx i dx j = 1,2,3). (15) 



At each point one can define a local scale factor a and a local Hubble parameter H 
by 

a 2 = (detTij)* ,H = - , (16) 

(X 

where dot denotes the time derivative. The Hubble parameter leads to the charac- 
teristic proper time on which the metric evolves. The characteristic comoving length 
on which it varies is denoted L : d^jk ~ L~Hj • The long wavelength approxi- 
mation is the assumption that the characteristic scale of spatial variation is much 
bigger than the Hubble radius, that is: 

-da jk < iij ^ aL > H- 1 . (17) 

(X 

Then we can drop the spatial curvature term in the Einstein equations in the low- 
est order. If we incorporate the curvature effect perturbatively, we get the under- 
standing of the non-linear evolution of the inhomogeneities. This direct method is 
technically useful. We presented the analysis in the Appendix for the purpose of 
the check. For the conceptual reason, here, we take another approach, i.e. , we 
consider the long-wavelength approximation in the context of the Hamilton- Jacobi 
formalism. 

Let us follow the method developed by Salopek and the co-orthors. They expand 
the generating functional in a series of terms according to the number of spatial 
gradients they contain: 

S = S^ + S^ + S^ + S^ + ---. (18) 

As a result the Hamilton- Jacobi equation can be solved perturbatively as we will 
show. The lowest order Hamilton- Jacobi equation is 

1 6S<®6SW t _. co + 1 _ _ , , 1 ,8S<-°\ 2 



■[likljl - , , , lijlkl\ + . I o\ rr (-rrY 



<P^1 Hi Hi J 2cu + 3 llJ '" J 2(2u + 3) ^ K 5<f> 

1 1 6S<®6SW j 5S^ 

2uj + 3 yff 5-fij 5<p v 5x 



It is difficult to obtain the complete solutions of the Hamilton- Jacobi equation. 
Fortunately, we need only the growing mode in the late stage or inflationary stage. 
Hence, we put the qusi-isotropic ansatz 

5(0) = _ 2 J dPxyfftHix) ■ (20) 

The diffeomorphyism invariance is automatically satisfied in this ansatz. Substitut- 
ing the ansatz (20) into eq.(19), we obtain 



Au + 6 dH 



H 2 = — — . (21) 

3u + 4 d X 

The solution of eq.(21) becomes 

where we ignored the integration constant, because it is a decaying mode. Hereafter, 
we choose the comoving gauge x,i — ■ In principle, we are free to choose the time 
slicing. For a dust fluid, a natural choice is a comoving gauge x,i — 0- Fortunately, 
in this case, we can also take the synchronous gauge, and x has a meaning of time. 
Equations of motion give 

* = 2^3* ' (23) 

* = Irl^- <24) 

Their solutions are 

= , (25) 

4^+4 

Hi = X 3 " +4 hij(x) , (26) 

where <f>(x) and hij(x) are arbitrary functions of spatial coordinates, which we call 
seed scalar and seed metric, respectively. 

The second order Hamilton- Jacobi equation is given by 

2u + 2 6SV> H 5S® 
2cu + 3 7ii + 2c; + 3^ 5<p 

+- y/i<t>R + (v^±Y J di<f>d j( f) + 2^A0 = . (27) 

ox <P 



To simplify the Hamilton- Jacobi equation, we utilize a conformal transformation of 
the three-metric and the Brans-Dicke scalar to define variables 



where Q and W satisfy 



fij = n~\j, (28) 
tj} = W' 1 ^ , (29) 



^ " + (30) 



dx 2lu + 3 
dW 1 
~dx 2uj + 3 



HW . (31) 



And then, 



„ 2u+2 

n = x —* , (32) 

W = x^ 1 ■ (33) 



Hence, the Hamilton- Jacobi equation reduces to 



5 X lfij '^ 
This is easily integrated to 



(34) 



S {2) = J d'x^j[4>R{l) - ^drfdrf] , (35) 

where we ignored the irelevant homogeneous solution. Now the equation of motion 



becomes 

= aJTa^f " 7 «S ] (36) 

H 1 11 ti 5SM 8S( 2 \ 
= 2^3 0+ 2^T3^ [0 ^~ 7 "'^7 ] ' ( } 

We expand <p as <p = 0^ + (f)^ H . 



Up to second order, the solution is given by 

2 ~ (3U + 4) 2 2^+6 . 

^ A ^ (5cu + 8)(2cj + 3)(2cj + 4) A K } 



where 



1 1 

F(/i, 0) = -0i?(/i) - JtwdiWrf + (2cu + 2)//''o : „ . (40) 
z z rt> 



Similarly, the metric is solved as 



^-^ V^V ^' (4i) 



where 



2w„7«7 2a; + 2, 1 , „- 2 ~ 

+ -^0^0 - 2u + 3 h ij^ h + ^<P;ij ■ 

The inhomogeneities of the space-time grow as ^- 2w + 4 / 2w + 3 i n this case. It is inter- 
esting to see if this tendency will continue or not. Therefore, we proceed to higher 
order calculations. Higher order generating functional is obtained by 

where 



R (2n) = £ 



n-l I SS (2k) SS (2n-2k) f u + 1 



hi 



likljl 2^ _|_ -j^Tij'Tfci 



(43) 



^1 1 ^ 6S (2k) 6S (2n-2k) n_l j j ^(2fc) ^(2n-2fc) 

+ ^ 2(2uj + 3)~/^^0 50 ^ 2^+3~/^ 7iJ ^~ 50 
Using the conformal transformation method, we obtain the recursion relation 

This leads to 



+ R {2n) = . (44) 



S {2n) = - j 'd 3 x £ dsxR™ [s X Mx) Jij(x)} 

3UJ + A { xj d 3 xR^ [ X , 0(x), 7y (x)] • (45) 



(2n + 3)uj + 4n + 4' 



From this expression, we can guess the following formal expansion 

^ = T,CnX— + —*\ (46) 

n 

where C n can be ,in principle, determined perturbatively. This formal expansion 
indicates the growing nature of the inhomogeneities of the space-time. It should be 
noted that this expression (46) coincides with that of the general relativity in the 
limit of uj — > oo. 

4 Cosmological Constant 

It is possible to show that the inhomogeneities grow or decay, as time increases, 
depending on the equation of the state for the perfect fluid matter. As we inves- 
tigated the dust matter which satisfies the strong energy condition, here, we will 
study the cosmological constant model as a typical one which does not satisfy the 
strong energy condition. The action is given by 

/ = J d 4 x^9[^R - ^d^dA - 2A] . (47) 

As is explained in the previous section, we can obtain the Hamilton- Jacobi equation 

1 5S 5S , u + l . 1 f 5S, 2 

-[liklji - , , , - lijlki] + , , ; , q \ /- (t") ( 48 ) 



ty/lHiSlu 2uj + 3 iij IKtl 2(2^ + 3)^/7 

' 1 ^'IT^TI - Vl^R + w^ytdiWrf + 2^A0 + 2A^ = . 



2w + 3 y /7" J 6'y ij 6<i> v ,r ' <p 
It is our task to solve the above Hamilton- Jacobi equation using the long-wavelength 
approximation. The lowest equation becomes 



1 ss^ssw 



UJ + 1 



+ 1 ,SS«>\ 2 



<t>y/l Hj hki L 2u> + 3" J H 2(2^ + 3)^7 5<P ' 

1 1 6SW6S<® 

z-Yij-z 7— + 2A^ = . (49) 



2^ + 3^7 Hj 

Here, we seek the quasi-isotropic solution again. For this purpose, we take the ansatz 



S(°> = -2J^/dV7tf(0- (50) 
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Substituting eq.(50) into eq.(49), we get 



0^ ~ 37? 90 ' ~ 90 + ~J ' ( ' 



Then, we obtain 



H=J^^. (52) 



Using the equations of motion, we obtain the lowest solution 

0(O) = (2c + 3xL + 5) t2 ' (53) 
if = t^h l3 {x). (54) 

If uj > 1/2, the space-time shows power-law expansion which is closely related to 
the extended inflationary scenario. The next order equation is 

. / nr ~ _i ss® r "1a - ~ ,5sw 

(2 ^ + 1) V(2c + 3)(6c + 5) 2 ^^~ + 2 V(2c + 3)(6^ + 5) 02 ^0" 

-v^0i? + ^7^09,0 + 2^A0 = . (55) 

The direct integration yields 
where 



4A 



V (2cj + 3)(6cj + 5) 
From the generating functional, the next order correction is calculated as 



(57) 



^ = - iu -i)^ + 7)Z+m^) ^" m - (58) 



7 « = "" ^ - 1)(2, + 3)(2. + 7) , '^ (ft) + (2, - 1)(2, + 7) ^ (,i)l - (59) 
In contrast to the dust case, the inhomogeneities will decay if uj > 1/2. This case 

corresponds to power law inflation. Recursive calculation gives the higher order 

solutions. It is not difficult to guess the expansion form for further correction: 

( f ) = J2t 2+n ^ 1 - 2 ^C n . (60) 

Therefore, the inhomogeneties of the gravitational constant will decay in the power 
law inflation case. 
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5 Discussions 



It is Dirac who first takes smallness of the following number seriously: 

Gm p m e ^ 40 
~ ru , 

e 2 

where G, e, m p and m e are the gravitational constant, the charge of the electron, the 
proton mass and the electron mass. Dirac' idea is simple; The gravitational constant 
varies with time like as 

Goc - . 

The smallness of the number merely indicates our Universe is old!! 

As a concrete model that realizes Dirac's idea, it is known that Brans-Dicke 
theory is most permissible in the sense that it appears as a low energy theory of 
superstring theory naturally and it passed the classical tests such as the equivalence 
principle. In the isotropic and homogeneous flat Universe, the time variation of the 
gravitational constant becomes 

G _ H{t) 1 

G ~ ~^TT K t ' 

where H is the Hubble parameter. 

Now we shall discuss the evolution of the gravitational constant in the inhomo- 
geneous case. First, let us consider the dust model and define the local Hubble, 

H = ^iHj 

2m + 2 1 {uo + 3)(3lu + 4) 2 ^ 

3uj + A X { 4(5cu + 8)(2cu + 3)(cu + i) X3 j ' 1 j 

where 

U = R(h) - Z-fPdM + -^\h kl 4>, kl . (62) 
<p z LO + 6 (p 

Then, within the accuracy of the present approximation, the rate of variation of 
gravitational constant is given by 
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G 
G 



H (3cu + 4) 2 ^±i^ (u 



LO + 1 



Q(5u + 8)(lu + 1)' 



where 



Q(h, 0) = - t-K^d^ + ^Ih^-M ■ 

(ft 2 (ft 



(63) 



(64) 



Due to the local curvature and the inhomogeneities of the Brans-Dicke field, the 
evolution of the effective coupling constant at that point is altered as in eq.(62). 
This effect becomes dominant in the late stage. For example, in the case = 0, 
the positive curvature enhanced the decreasing effect, and the negative curvature 
surpressed the decreasing effect. In general, the inhomogeneities of the gravitational 
constant itself gives influence to the evolution rate of G. Thus the inhomogeneities 
becomes significant. Let us look at the evolution of the density of the dust 



_ i 
7 1 



1 6„+6 r (cj + 3)(3cj + 4) 2 
= y"^+4{l -I 1 il 1 

Vh 1 4(5^ + 8)(2cu + 3)(u + 2) 



2^+4 



(65) 



This shows that high density region and the small gravitational constant region 
coincides. This may give some interesting consequences to astrophysics. 

In the case of the cosmological constant, the time evolution of the gravitational 
constant is given by 



G 
G 



2uj + 1" 



-H 



(Au- l)(4u; + 3) 



6(u - l)(2w + l)(2w + 3)(2cu + 7) 



l-2u; 



R(h) 



(66) 



where 



H = 



2lu + 1 
2t 



2u(u + 2) 



-t 



1-2U! 



R(h) 



(67) 



3(lu - l)(2w + l)(2u + 3)(2cj + 7) 

In the case of the power law inflation, the homogenization of the time variation rate 
of gravitational constant occurs. 
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6 Conclusion 



As a fundamentally important theory, we have studied the Brans-Dicke theory using 
the long- wavelength approximation. First of all, we presented the Hamilton- Jacobi 
equation for the Brans-Dicke theory which has different structure from the that 
of general relativity. The Hamilton- Jacobi equation thus obtained is considered as 
the semi-classical equation correspondent to the Wheeler-DeWitt equation for the 
Brans-Dicke theory. Applying the method of Salopek and co-orthers, we have ob- 
tained the approximate solutions up to the first order of curvature correction in the 
case of dust fluid and the cosmological constant. As the Brans-Dicke scalar can 
be regarded as the effective gravitational constant, the non-linear evolution of the 
inhomogeneities of the Brans-Dicke field is interesting. In the above two cases, we 
investigated this problem. From the results, we can conclude that the inhomogeni- 
eties of the gravitational constant will decay in the case of inflatinary matter and 
grow in the case of ordinary matter which satisfies the dominant energy condition. 
We have also calculated the approximate solutions using the direct method of Comer 
et.al. in the Appendix. 
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APPENDIX: Anti-Newtonian Formalism 

The equations of motion for the Brans-Dicke theory is 



G 



Tmhv + -t(4> 




9^) + j 2 (d^d v <\> - -g^d^d^) (68) 

(69) 



□ = 



14 



where Tm^v is the energy momentum tensor of the matter, hereafter we will con- 
sider the irrotational dust fluid. In the synchronous gauge, we can write the above 
equations in the following way: 



d_ 
dt 



K + K\K\ 



-KL + K\i 



d 



KK{ + -Kl +< 3 ) R{ 



-0 - K<p + A0 



1 +i^ )fe -^0 



i 



uo 

-2' 



47T 



1 





1 



0£-i<10 + ^(-0-^ + A0) 



+-d*<f>d k <f> + — P 5\ 



Air 
J 



8tt 



2u + 3 



P 



(70) 
(71) 

(72) 
(73) 



where p is the energy density of the dust fluid. 
Define 

K] = Y? 3+ l -5)K 
We can write the above equations in the following way: 



d_ 

dt 



K + -K 2 



i^tf + l R 



1(3) 



3 u 



uo 1 



<9t 3 



= --0 2 -(a; + 3)^ 



(74) 



(75) 



(76) 



2 " 

-(cu + 2)K^ + (cu + 2)iA0 (77) 



+^(5>d fc 0-I^ m 0d m 0) (78) 

We now consider these equations order by order in the gradient expansion. At 
zeroth order, we neglect the terms which have twice derivative, i.e., ^R,A(f) etc.. 
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Futhermore, we impose quasi-isotropic nature E*- = 0. In this order, we obtain the 
solutions 



4^+4 

~ 2 
= 0(x)t 3 "+ 4 , 

p = 



x) Auj + 6 e^+e 

t 3^+4 . 



1 

3 



<5]A0) 



(79) 
(80) 
(81) 



8tt 3. + 4 

Here, /iy and depend on the spatial coordinates. Substituting (79) ~ (81) into 
eqs.(75)~ (78) and keping the next oeder terms, we obtain 



(82) 



(3. + Af 



(5. + 8)(2. + 4)(2. + 3) 



2^+6 



-1 

-<PR + 2(u+l)A<f ) -^-d k <f ) d k < 
2 2<p 



(83) 



K 



i(2) 



3a; + 4 
5c7+8 



y 2(2. + 3) J 
. + 1 



1 ^ 



2. + 3 J 



a; 



2(2. + 3)^^ 



(84) 



(5. + 8)(. + 2) L J 2(2. + 3) 



2. + 3 



1 :<^W] 



JV 2(2. + 3p0 2 ' 
The above results completely agree with the Hamilton- Jacobi calculation. 



(85) 
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